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Abstract 



The prediction of the strong coupling assuming (supersymmetric) coupling 
constant unification is reexamined. We find, using the new electroweak data, 
a s (Mz) ~ 0.129 ±0.010. The implications of the large a s value are discussed. 
The role played by the Z beauty width is stressed. It is also emphasized that 
high-energy (but not low-energy) corrections could significantly diminish the 
prediction. However, unless higher-dimension operators are assumed to be 
suppressed, at present one cannot place strong constraints on the super-heavy 
spectrum. Non-leading electroweak threshold corrections are also discussed. 
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Assuming the minimal supersymmetric extension of the standard model (MSSM) |J 
between the weak and some high scale, one finds 0] that the extrapolated electroweak and 
strong couplings approximately unify at a scale M G ~ 3 x 10 16 GeV (the grand unification 
scale). Alternatively, assuming coupling constant unification, one can use the precisely 
measured weak angle s 2 (Mz) and fine-structure constant a(Mz) to predict the Z-pole strong 
coupling a s (Mz). Model-dependent corrections are typically of order 10%, i.e., comparable 
to the experimental uncertainty in a s (Mz), and need to be included consistently H. Below, 
we update and extend our discussion of the a s (M z ) prediction H3H6]. We find that for the t- 



quark pole mass m?° le > 160 GeV, the positive corrections proportional to m 2 are sufficiently 
large that the sum of the (Yukawa, threshold, and operator) model-dependent corrections 
must cancel or be negative for unification to hold. Ignoring possible high-scale matching 
corrections, tan /3 « 1 and heavy superpartners are preferred (tan (3 = (H 2 ) / (Hi)). However, 
large negative high-scale threshold and nonrenormalizable operator (NRO) corrections are 
possible. The former depend on the details of the grand- unified theory (GUT), while the 
latter |7j are gravitationally induced and are generic. Below, we review our formalism and 
discuss our results and their implications. We also comment on non-logarithmic superpartner 
corrections, implications of the anomalous Z — ► bb width, extended models, and on various 
aspects of the large QCD coupling. A comprehensive analysis is presented in Ref. M. 
The prediction for a s (M z ) readsQ 



a s {M z ) = a° L (M z ) + 0.014 + H as + 
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where s 2 (M z ) is the true (MS) weak angle and Sq is the value it would have for mf 
The 0.014 correction is a (model-independent) two-loop gauge correction and the function 
H ag is a smaller (model-dependent) two-loop Yukawa correction. a 2 /28n is a finite scheme- 
dependent term. The model-dependent function A as sums threshold and NRO corrections at 
low and high scales. Substituting in (0) the (MS) input values a(M z ) = 1/(127.9±0.1) 

and 



1 Hypercharge is properly normalized, i.e., s 2 (Mg) = 3/8. 

2 We do not explicitly treat smaller logarithmic dependences on rr?° le . They are included in the 
uncertainty. The 0.88 factor incorporates higher-order QCD corrections which were not included 
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s 2 (M z ) = 0.2316 ± 0.0003, (4a) 

rrC = 160^ + 13 In ^ (4b) 

(rriho is SM-like light Higgs boson mass0), one has (in the MS scheme) 

a s (M z ) - A as = 0.129 ± 0.001 + 3.1 x 10~ 7 GeV~ 2 [(m p t ole ) 2 - (160 GeV) 2 ] + H Qs . (5) 
The higher values of mF° le (e.g., compared to ||) and lower value of the weak angle implied 



by recent data ||TT| increase the predicted a s . An even higher central a s value of 0.130 would 
be predicted for the value mF t ole = 174 ± 16 GeV suggested by the CDF t-quark candidate 
events Two- loop Yukawa corrections are negative but are typically negligible. They 

can be important if the Yukawa couplings of the t and/or 6-quark, h t and h^, respectively, 
are large, i.e., for tan/5 ~ 1 or tan/? > 50. We find |5| 

- 0.003 < H as (h t , h) = H as (m p t ole Mn(3) < 0. (6) 

For h t ~ max[h t (m t )] ~ 1.1 (and h b ~ 0) one has || H as ~ —0.1 xa 2 xh 2 ~ —0.002. In gen- 
eral, one can substitute a one-loop semi-analytic expression for h 2 and integrate iteratively 
|i~3| (a similar procedure leads to our result for the gauge two- loop correction @.|[). 

The coupling constant unification is shown in detail in Fig. 1 for various values of 

a s (M z ) = 0.12 ± 0.01 and for A Qs = and H as 0.0005. In the absence of threshold 

corrections, and for reasonable mf°' e , coupling unification requires a s (M z ) ^ 0.127. Below, 
we show that typically |A a J < 0.01. Thus, we obtain from coupling constant unification, as- 
suming no conspiracies among different model-dependent corrections, a s (Mz) ^ 0.12. This 
is in a good agreement with Z-pole extractions of a s , but is slightly higher than some ex- 
tractions based on deep inelastic scattering (DIS) and quarkonium spectra. The prediction 
is compared with the data in Table 1 (from (T^j). The a s measurement and the possibility of 



light gluinos (that correct the a s extrapolation between the quarkonium and weak scales by 
~ 10%) are further discussed in Ref. |15],[l6|. We note, in passing, that light colored scalars 
would correct the a s extrapolation negligibly, i.e., a light scalar top would affect the extrap- 
olation of a s measured at low-energy to the Z-pole by less than 1%. Models (in particular, 
NRO's) can be constructed with large (> 10 — 20%) and negative GUT scale contributions 
to A as . Such models would violate our no-conspiracy assumption, but cannot be excluded. 
Hence, even if supersymmetry is characterized by experiment and the superpartner contri- 
bution to A aB (see below) is found to be positive, coupling constant unification will not 



3 The authors of perform a best fit to all W, Z and neutral current data assuming 60 < 
nifto < 150 GeV with a central value m^o = Mz for the SM-like light Higgs boson mass. (Other 
possible light particle corrections are discussed separately below.) In the (non-supersymmetric) 
standard model one assumes a larger Higgs mass range 60 < m h o < 1000 GeV with a central value 
of 300 GeV. This leads to the prediction m^° le = 175 ± ll^ig GeV, where the second uncertainty 
is from m h o. 
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be completely ruled out even for a s (Mz) ~ 0.11. However, one will be able to sufficiently 
constrain GUT's only if the superpartner contribution is large and positive (i.e., if NRO's 
with perturbative coefficients are not sufficient to rectify the prediction). 

The situation in the non-supersymmetric extension is quite different since (a) su- 
persymmetry doubles the GUT sector, (b) NRO's are typically suppressed in the non- 
supersymmetric case by powers of (Mq / 'Mpi anc k) ~ 10~ 5 , and (c) the corrections oc ct 2 s (Mz) 
are suppressed by a ~ (0.07/0.13) 2 factor in comparison to the MSSM One can rectify 
this situation by considering large logarithms and/or certain complicated chain-breaking 
scenarios with additional particles, i.e., intermediate scales (which, however, could be con- 
structed to be C(10 16 GeV) [[FT] or 0(1 TeV) The predictive power of a desert theory 
is lost in such a case. 

Next, we discuss in greater detail the possible model-dependent contributions to the 
(9(10%) correction function 

— 19a 2 Msusy 

A as & — — — - In — — h GUT threshold corrections + NRO corrections. (7) 

287r Mz 

The parameter Msusy 0] is a weighted sum of all superpartner and heavy Higgs boson 
mass logarithms which determines the (leading-logarithm) contribution to A ag || [A Qs ~ 
-0.003 \n(M S usY/M z )]. Specifically, 

Msusy = U^ 3s[ " 2 " - (8) 

i 

where the index i runs over all superpartner and heavy Higgs particles. We defined the 
/?-function coefficients b % - = a(Si)Ni(j)ti(j), where a(Si) = |, |, — y for a particle i of spin 
Si = 0,1,1, respectively, Ni(j) is the appropriate multiplicity, and t;(j) = 0, |, 2, 3, |(^) 2 
for a singlet, a particle in the fundamental representation of SU(N), an SU(2) triplet, an 
SU (3) octet, and for j = 1 and a particle with hypercharge Y, respectively. Because of mass 
non-degeneracies between colored particles (whose masses are sensitive to the gluino mass), 
the Higgs and Higgsino particles (whose masses are sensitive to /i), and the scalar leptons 
(whose masses are sensitive to scalar mass boundary condition), and because of the different 
weights assigned to the different particles, Msusy is not simply the geometric mean of the 
rrii. In particular, the negative powers in (§) imply that Msusy can be (and generally is) 
much smaller than the actual masses of the superpartners. In Fig. 2 we calculate Msusy 
for more than a thousand arbitrary^ MSSM's which are consistent with the electroweak 
symmetry breaking, a neutral lightest supersymmetric particle, and sparticle masses above 
experimental lower bounds and below ~ 2 TeV (see [[TJ,^(J). M S usy is proportional to the 
Higgsino mass parameter fi f2l[ and is indeed lower than the actual superpartner and Higgs 



boson masses. From Fig. 2 we have the approximate upper bound Msusy ~ 250 — 300 GeV 
(or the lower bound A^ SY > -0.003). 

As mentioned above, H ag is large and negative for tan/3 ~ 1. Also, M S usy oc \fi\ oc 
l/[tan 2 /3 — 1] is maximized in that region of the parameter space {Msusy is shown as 



4 We assume universality of the soft parameters at Mq. However, similar results for Msusy hold 
in more general scenarios. 
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a function of tan/3 in Fig. 3). The proportionality factor depends on and grows with the 
superpartner masses. Thus, a heavy spectrum and tan /3 ~ 1 are slightly preferred. This 
observation is consistent with b — r Yukawa unification (which we do not require here), which 
is constrained by the interplay between the large predicted values of a s and the Yukawa- 
unification preference of moderate a s values ||. (The large QCD radiative corrections to hf, 
constrain one to regions of the parameter space in which large Yukawa coupling can partially 
compensate for these corrections^].) In that region one has the spectacular constraint on the 
Higgs boson mass m h o < 100 (110) GeV for mF t ole < 160 (175) GeV at one loop (and a 
stronger bound applies at two loops) 



It was recently suggested that the Z-pole couplings should be extracted from the data 
assuming the full MSSM p4|| . This is the case if the model contains some particles (aside 
from the SM-like Higgs boson) lighter than ~ 100 — 150 GeV. However, assuming the heavy 
MSSM limit, SU(2) breaking mixing and other non-leading effects are negligible and our 
leading-logarithm formula, which is derived using renormalization-group techniques, is an 
excellent approximation. Otherwise, light particle (non-logarithmic) effects can be accounted 
for in the same manner used to describe the quadratic m t dependence ||[)|||, i.e., by the 
perturbative expansion^ 



s\M z ) = 4{M Z ) + u ; - ArT + Arf SY , (9) 

where Ar z sums (universal) corrections to the Z-boson mass M z . The leading contri- 
bution to Ar z p is given in eq. (^) and Ar| USY has been calculated in Ref. [26|. In fact, it 



„SUSY 



is useful to subtract from Ar| USY leading logarithms summed by Msusy and reserve Ar| u 
to denote only additional contributions of light superpartners. The correction function (|7|) 
is modified accordingly, A ag — > A ag — 1.16Ar| USY . The different contributions to Ar| USY 
are correlated in a given model, and their interplay determines its magnitude and overall 
sign. We find p7| that non- logarithmic corrections typically conspire with the term and 



increase the a s prediction, in some cases, by a few percent. Thus, heavy superpartners are 
preferred beyond the leading order. 

On a similar note, it has been observed that if supersymmetry significantly modifies 
the Z hadronic width (so that the Z —>■ bb anomaly is accounted for) then a s extracted 
from the Z line shape is diminished significantly (e.g., 0.126 — > 0.112) UTTf , and this effect 
was even promoted as a possible resolution of the discrepancy between low and high-energy 
extractions of a s |16| . Such a scenario would require either light Higgsinos and large Yukawa 
couplings or a very large tan/3 and a light pseudo-scalar Higgs boson p8| , i.e., |/i| ^ 0(M Z ). 



However, a scheme with a small fx parameter is not favored in GUT models p0| . From 
our discussion above it is also clear that a solution involving light Higgsinos (or a light 
pseudo-scalar) is strongly disfavored by the a s prediction: 



5 Finite superpartner loops [^] modify only the allowed large tan/3 region. 

6 One could calculate the corrections to all fitted observables, or risk a minor inconsistency and 
calculate only (universal) corrections to the input parameter (M z in our case). The latter scheme, 
which we follow, is sufficient for our current purposes. 
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1. The extracted a s line-shape value would decrease (in agreement, however, with low- 
energy extractions). 

2. The predicted ot s value would increase due to leading-logarithm [oc — ln(|/j|/Mz)] and 
possibly non-logarithmic threshold corrections. 

3. The central value of the fitted m? t ole [eq. ([£51)1 would grow to 163 GeV, further increas- 
ing the a s prediction by 0.0003. 

Thus, the Z — > bb anomaly, if not resolved, contains strong implications for supersymmetric 
models and could even rule out the simplest and most attractive unification scenarios. 

Lastly, we consider possible high-scale contributions to the correction function. Unlike 
the MSSM case, in which the particles and their mass range are dictated by the model, 
the details of the high-scale corrections are ambiguous. In the minimal SU(5) model |29| 



negative threshold corrections in (|7|) due to super-heavy color triplet Higgs supermultiplets 
are strongly constrained by the non-observation of proton decay |j30"| , and the GUT-scale 
threshold correction contribution to A as is typically positive. (This observation, however, 
need not hold in extended models.) Nevertheless, one cannot extract strong constraints on 
the GUT spectrum. Gravitationally induced operators (suppressed by M G /M P i anck ~ 0.001) 
split the M G gauge couplings (in a correlated manner) and correct the a s prediction in 
proportion to their effective strength, rj, which is a free parameter and can have either sign. 
One has[] A^ RO pa 0.00 57?. Constraining the NRO corrections to stay perturbative so that 
the calculation is consistent (higher-order terms are negligible in this case) one has \r)\ ^ 2 
(\r)\ ~ 3 is an extreme but still acceptable choice). Thus, NRO's with a non-negligible 
and negative r\ could smear light and heavy threshold corrections. Unless r\ ^ and/or 
Msusy Mz (which could also imply positive non- logarithmic corrections), no significant 
constraints can be placed on the super-heavy spectrum at present. On the other hand, the 
minimal SU(5) model (where threshold corrections are strongly constrained) would require 
NRO's with T) < if a s (Mz) ^ 0.125. (A similar observation was made recently in Ref. 
PHI .) Thus, unification and quantum gravity may be inseparable. 

Regarding the unification scale, corrections that increase the unification scale would 
typically also increase the prediction for a s {Mz) p2]| , and are thus difficult to construct 
[in particular, for s 2 (Mz) ~ 0.2316]. This is true for contributions to A aB as well as for 
an additional matter family \a s (Mz) — > 0.132] or additional pairs of Higgs doublets [which 
lead to non-perturbative values of a s (M z )]- This is easily understood if we write a° L as a 
function of the unification scale M and of a(M z ) (see Fig. 4), 



o?{M z ) = MMZJ ; , (10) 
s V ' 3-60a(M z )t' v ; 



where t = (l/2vr) ln(M/M z ). Naively substituting, e.g., M = M string ~ 5 x 10 17 GeV {53 



one has a s (Mz) > 0.2. By carefully adjusting operator and super-heavy threshold correction 



7 The proportionality factor is calculated here in the SU(5) theory [fj], and its normalization is 
different by a factor of four than in |3| . 
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contributions to A as , one could increase Mq by an order of magnitude while maintaining an 
acceptable prediction for a s [f|||. However, in general, to rectify the string and unification 
scales (in level-one models) one has to compromise the predictive power of the unification 
scenario 



34 so that the correlation between a s and t is modified. 



To conclude, we have shown that typically one expects a large QCD coupling in super- 
symmetric unified models (and even more so when considering a typical MSSM spectrum). 
This constitutes an interesting signature and has implications for, e.g., Yukawa unification, 
corrections to the unification scale, and the overall sign of the correction function A ag , and 
is in possible conflict with low-energy data. However, it does not yet allow a significant 
constraint on the super-heavy spectrum because of possible gravitational corrections. We 
also pointed out the interesting role that the Z hadronic width might play in supersymmet- 
ric GUT's, and suggested a simple formula that extends our treatment of mf°' e -dependent 
electroweak corrections to the supersymmetric sector. 
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TABLES 



TABLE I. Values of a s (Mz) extracted from different processes (and extrapolated to Mz if 
relevant). The different values are ordered according to the energy scale of the relevant process. 

Bjorken sum rules 0.122+°;°^ 

r -» hadrons (CLEO) 0.114 ± 0.003 

r -» hadrons (LEP) 0.122 ± 0.005 

Deep inelastic scattering 0.112 ± 0.005 

J/* (lattice) 0.110 ±0.006 

T (lattice) 0.115 ±0.002 

T, J/* (decays) 0.108 ± 0.010 

ep -> 2 + 1 jet rate (HERA) 0.121 ±0.015 

e+e" event shape (LEP) 0.123 ± 0.006 

Z line shape (LEP) 0.126 ± 0.005 
Prediction 0.13 ±0.01 
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FIGURES 



FIG. 1. MSSM evolution of a± i 2 (solid lines) and of 03 (dashed lines) in the vicinity of the a>i : 2 
unification point (the scale M is in GeV). a s (M z ) = 0.110, 0.115, 0.120, 0.125, 0.130; m p t ole = 160 
GeV; tan (3 = 4; and A as = 0. 

FIG. 2. Msusy as a function of the /i parameter. The different universal soft parameters 
and tan/3 are picked at random in the allowed parameter space (see text). m pole = 160 GeV. 
Msusy = Mz is denoted for comparison. (All masses are in GeV.) 



FIG. 3. Same as in Fig. 2 except a function of tan/3. 



FIG. 4. The Z-pole weak angle and strong coupling are predicted as a function of the unifica- 
tion scale M. A given value of s 2 (Mz) corresponds to a fixed choice for M, e.g., s 2 (Mz) = 0.2359 
corresponds to M = 10 16 GeV. MSSM /3-functions are assumed. Two-loop Yukawa corrections are 
taken into account assuming m p t ole = 160 GeV and tan/3 = 4. (A Qs = 0.) s 2 (M z ) = 0.2316±0.0003 
and a s (Mz) = 0.12 ± 0.01 are indicated for comparison. 
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Assuming the minimal supersymmetric extension of the standard model (MSSM) [1] 
between the weak and some high scale, one finds [2] that the extrapolated electroweak and 
strong couplings approximately unify at a scale Mq ~ 3 X 10 16 GeV (the grand unification 
scale). Alternatively, assuming coupling constant unification, one can use the precisely 
measured weak angle s 2 (Mz) and fine-structure constant a(Mz) to predict the Z-pole strong 
coupling a s (Mz)- Model-dependent corrections are typically of order 10%, i.e., comparable 
to the experimental uncertainty in a s (Mz) } and need to be included consistently [3]. Below, 
we update and extend our discussion of the a s (Mz) prediction [3-6]. We find that for the t- 
quark pole mass m p t ole > 160 GeV, the positive corrections proportional to m 2 are sufficiently 
large that the sum of the (Yukawa, threshold, and operator) model-dependent corrections 
must cancel or be negative for unification to hold. Ignoring possible high-scale matching 
corrections, tan f3 ~ 1 and heavy superpartners are preferred (tan f3 = (H 2 ) / (Hi)). However, 
large negative high-scale threshold and nonrenormalizable operator (NRO) corrections are 
possible. The former depend on the details of the grand-unified theory (GUT), while the 
latter [7] are gravitationally induced and are generic. Below, we review our formalism and 
discuss our results and their implications. We also comment on non-logarithmic superpartner 
corrections, implications of the anomalous Z — >■ bb width, extended models, and on various 
aspects of the large QCD coupling. A comprehensive analysis is presented in Ref. [8]. 

The prediction for a s (Mz) reads 1 

a s (Mz) = af{M z ) + 0.014 + H a . + + 3.1 x 10" 7 GeV" 2 \(mf e f - (mff) 2 ] + A Q 

Z07T L J 



where 



a^(Mz) 



7a(M : 



z 



15s§(M z )-3' v ; 

is the lowest order prediction, and 2 

s 2 {Mz) = s 2 (M z ) - 0.88 x 10" 7 GeV" 2 \{m v t ole f - {m v t ° le f \ , (3) 



pole pole 



m f = m 



to 



where s 2 (Mz) is the true (MS) weak angle and Sq is the value it would have for 
The 0.014 correction is a (model-independent) two-loop gauge correction and the function 
H as is a smaller (model-dependent) two-loop Yukawa correction. a 2 /287r is a finite scheme- 
dependent term. The model-dependent function A„, sums threshold and NRO corrections at 
low and high scales. Substituting in (1) the (MS) input values [9-11] a(Mz) = 1/(127. 9±0.1) 
and 



1 Hypercharge is properly normalized, i.e., s 2 (Mq) = 3/8. 

2 We do not explicitly treat smaller logarithmic dependences on m^ ole . They are included in the 
uncertainty. The 0.88 factor incorporates higher-order QCD corrections which were not included 
in [3]. 
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s 2 (Mz) = 0.2316 ± 0.0003, (4a) 
m-' e = 160^ + 131n^ (4b) 

(m h o is SM-like light Higgs boson mass 3 ), one has (in the MS scheme) 
a s (M z ) - A a , = 0.129 ± 0.001 + 3.1 x 10" 7 GeV" 2 \{m pole ) 2 - (160 GeV) 2 ] + H a ,. (5) 

The higher values of m p t ole (e.g., compared to [3]) and lower value of the weak angle implied 
by recent data [11] increase the predicted a s . An even higher central a s value of 0.130 would 
be predicted for the value m pole = 174 ± 16 GeV suggested by the CDF t-quark candidate 
events [12]. Two-loop Yukawa corrections are negative but are typically negligible. They 
can be important if the Yukawa couplings of the t and/or 6-quark, h t and respectively, 
are large, i.e., for tan f3 ~ 1 or tan f3 <; 50. We find [5] 

-0.003 < H a ,(h t , h) = H as {mf e Mnfi) < 0. (6) 

For h t ~ max[/i((m ( )] ~ 1.1 (and h b ~ 0) one has [3] H as ~ —0.1 X a 2 s X h 2 ~ —0.002. In gen- 
eral, one can substitute a one-loop semi-analytic expression for h 2 and integrate iteratively 
[13] (a similar procedure leads to our result for the gauge two-loop correction [3,8]). 

The coupling constant unification is shown in detail in Fig. 1 for various values of 

a s (M z ) = 0.12 ± 0.01 and for A„, = and H a , 0.0005. In the absence of threshold 

corrections, and for reasonable m pole , coupling unification requires a s (Mz) ^ 0.127. Below, 
we show that typically |A„J < 0.01. Thus, we obtain from coupling constant unification, as- 
suming no conspiracies among different model-dependent corrections, a s (Mz) ^ 0.12. This 
is in a good agreement with Z-pole extractions of a s , but is slightly higher than some ex- 
tractions based on deep inelastic scattering (DIS) and quarkonium spectra. The prediction 
is compared with the data in Table 1 (from [14]). The a s measurement and the possibility of 
light gluinos (that correct the a s extrapolation between the quarkonium and weak scales by 
~ 10%) are further discussed in Ref. [15,16]. We note, in passing, that light colored scalars 
would correct the a s extrapolation negligibly, i.e., a light scalar top would affect the extrap- 
olation of a s measured at low-energy to the Z-pole by less than 1%. Models (in particular, 
NRO's) can be constructed with large (> 10 — 20%) and negative GUT scale contributions 
to A„ s . Such models would violate our no-conspiracy assumption, but cannot be excluded. 
Hence, even if supersymmetry is characterized by experiment and the superpartner contri- 
bution to A Qs (see below) is found to be positive, coupling constant unification will not 
be completely ruled out even for a s (Mz) ~ 0.11. However, one will be able to sufficiently 
constrain GUT's only if the superpartner contribution is large and positive (i.e., if NRO's 
with perturbative coefficients are not sufficient to rectify the prediction). 



3 The authors of [11] perform a best fit to all W, Z and neutral current data assuming 60 < m h a < 
150 GeV with a central value m h a = Mz for the SM-like light Higgs boson mass. (Other possible 
light particle corrections are discussed separately below.) In the (non-supersymmetric) standard 
model one assumes a larger Higgs mass range 60 < m h o < 1000 GeV with a central value of 300 
GeV. This leads to the prediction m P ° 6 = 175 ± H^Jg GeV, where the second uncertainty is from 
m h o. 
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FIG. 1. MSSM evolution of a.\^ (solid lines) and of 013 (dashed lines) in the vicinity of the 



unification point (the scale M is in GeV). a s (M z ) = 0.110, 0.115, 0.120, 0.125, 0.130; m 
GeV; tan/3 = 4; and A„, = 0. 



pole 



160 




The situation in the non-supersymmetric extension is quite different since (a) su- 
persymmetry doubles the GUT sector, (b) NRO's are typically suppressed in the non- 
supersymmetric case by powers of (Mq / Mpi anc k) ~ 10~ 5 , and (c) the corrections oc a 2 s (Mz) 
are suppressed by a ~ (0.07/0. 13) 2 factor in comparison to the MSSM [3,8]. One can rectify 
this situation by considering large logarithms and/or certain complicated chain-breahing 
scenarios with additional particles, i.e., intermediate scales (which, however, could be con- 
structed to be (9(10 16 GeV) [17] or Oil TeV) [18]). The predictive power of a desert theory 
is lost in such a case. 

Next, we discuss in greater detail the possible model-dependent contributions to the 
0(10%) correction function 

-19a 2 Msusy 

A„ - In h GUT threshold corrections + NRO corrections. (7) 

28tt M z J 

The parameter Msusy [3] is a weighted sum of all superpartner and heavy Higgs boson 
mass logarithms which determines the (leading-logarithm) contribution to A„, [3] [A„, ~ 
-0.003 ln(M S usY/M z )]. Specifically, 

- — \Ab l -— b' + — b'] 

Msusy = Iim t ss[ 1 + (8) 



where the index i runs over all superpartner and heavy Higgs particles. We defined the 

b) = a(S t )N t (j)t t (j) } where a(S t ) = |,|,-^ 



/3-function coefficients b % - = a(Si)Ni(j)ti(j), where a(Si) = ^, |, — ^ for a particle i of spin 



Si = 0,^,1, respectively, Ni(j) is the appropriate multiplicity, and ti(j) = 0,|,2,3,|(y) 2 
for a singlet, a particle in the fundamental representation of SU(N) } an SU{2) triplet, an 
SU{3) octet, and for j ' = 1 and a particle with hypercharge Y, respectively. Because of mass 
non-degeneracies between colored particles (whose masses are sensitive to the gluino mass), 
the Higgs and Higgsino particles (whose masses are sensitive to //), and the scalar leptons 
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TABLE I. Values of a s (Mz) extracted from different processes (and extrapolated to Mz if 
relevant). The different values are ordered according to the energy scale of the relevant process. 



Bjorken sum rules 


n i 99 +0.005 

U.±^_Q Q 09 


t -> hadrons (CLEO) 


0.114 ± 0.003 


t — ► hadrons (LEP) 


0.122 ±0.005 


Deep inelastic scattering 


0.112 ± 0.005 


«//$ (lattice) 


0.110 ±0.006 


T (lattice) 


0.115 ±0.002 


T, J/$ (decays) 


0.108 ±0.010 


ep -+ 2 + f jet rate (HERA) 


0.121 ± 0.015 


e + e~ event shape (LEP) 


0.123 ±0.006 


Z line shape (LEP) 


0.126 ±0.005 


Prediction 


0.13 ±0.01 



(whose masses are sensitive to scalar mass boundary condition), and because of the different 
weights assigned to the different particles, Msusy is n °t simply the geometric mean of the 
rrii. In particular, the negative powers in (8) imply that Msusy can t> e (and generally is) 
much smaller than the actual masses of the superpartners. In Fig. 2 we calculate Msusy 
for more than a thousand arbitrary 4 MSSM's which are consistent with the electroweak 
symmetry breaking, a neutral lightest supersymmetric particle, and sparticle masses above 
experimental lower bounds and below ~ 2 TeV (see [19,20]). Msusy is proportional to the 
Higgsino mass parameter // [21] and is indeed lower than the actual superpartner and Higgs 
boson masses. From Fig. 2 we have the approximate upper bound Msusy ~ 250 — 300 GeV 
(or the lower bound A^ SY > -0.003). 

As mentioned above, H as is large and negative for tan/? ~ 1. Also, Msusy ck \fi\ oc 
^l/[tan 2 {3 — 1] is maximized in that region of the parameter space (Msusy is shown as 
a function of tan {3 in Fig. 3). The proportionality factor depends on and grows with the 
superpartner masses. Thus, a heavy spectrum and tan {3 ~ 1 are slightly preferred. This 
observation is consistent with b — r Yukawa unification (which we do not require here), which 
is constrained by the interplay between the large predicted values of a s and the Yukawa- 
unification preference of moderate a s values [8]. (The large QCD radiative corrections to hj, 
constrain one to regions of the parameter space in which large Yukawa coupling can partially 
compensate for these corrections 5 .) In that region one has the spectacular constraint on the 
Higgs boson mass ra^o < 100 (110) GeV for m p t ole < 160 (175) GeV at one loop (and a 
stronger bound applies at two loops) [23,19,8]. 

It was recently suggested that the Z-pole couplings should be extracted from the data 
assuming the full MSSM [24]. This is the case if the model contains some particles (aside 
from the SM-like Higgs boson) lighter than ~ 100 — 150 GeV. However, assuming the heavy 



4 We assume universality of the soft parameters at Mq. However, similar results for Msusy hold 
in more general scenarios. 

5 Finite superpartner loops [22] modify only the allowed large tan (3 region. 



5 



51 

51 — — — 11 



51 



51 

51 - 




-2000 -7500 -7000 -500 500 7000 7500 2000 

FIG. 2. M5[/5y as a function of the fj, parameter. The different universal soft parameters 
and tan/? are picked at random in the allowed parameter space (see text). e = 160 GeV. 
Msusy = Mz is denoted for comparison. (All masses are in GeV.) 

MSSM limit, SU{2) breaking mixing and other non-leading effects are negligible and our 
leading-logarithm formula, which is derived using renormalization-group techniques, is an 
excellent approximation. Otherwise, light particle (non-logarithmic) effects can be accounted 
for in the same manner used to describe the quadratic m t dependence [3,6,8], i.e., by the 
perturbative expansion 6 

s*(m z ) = 4(m z ) + s \ {1 ~f l Ar T + A 4 USY ] , (9) 

where Ar z [25] sums (universal) corrections to the Z-boson mass M Z - The leading contri- 
bution to Ar^ p is given in eq. (3) and Ar|? SY has been calculated in Ref. [26]. In fact, it 
is useful to subtract from Ar|? SY leading logarithms summed by Msusy an d reserve Ar|? SY 
to denote only additional contributions of light superpartners. The correction function (7) 



One could calculate the corrections to all fitted observables, or risk a minor inconsistency and 
calculate only (universal) corrections to the input parameter (M z in our case). The latter scheme, 
which we follow, is sufficient for our current purposes. 
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FIG. 3. Same as in Fig. 2 except a function of tan/?. 
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is modified accordingly, A„, — > A„ s — 1.16Ar^ USY . The different contributions to Ar|? SY 
are correlated in a given model, and their interplay determines its magnitude and overall 
sign. We find [27] that non-logarithmic corrections typically conspire with the m\ term and 
increase the a s prediction, in some cases, by a few percent. Thus, heavy superpartners are 
preferred beyond the leading order. 

On a similar note, it has been observed that if supersymmetry significantly modifies 
the Z hadronic width (so that the Z — s- bb anomaly is accounted for) then a s extracted 
from the Z line shape is diminished significantly (e.g., 0.126 — > 0.112) [11], and this effect 
was even promoted as a possible resolution of the discrepancy between low and high-energy 
extractions of a s [16]. Such a scenario would require either light Higgsinos and large Yukawa 
couplings or a very large tan f3 and a light pseudo-scalar Higgs boson [28], i.e., \fi\ < 0(Mz). 
However, a scheme with a small // parameter is not favored in GUT models [20]. From 
our discussion above it is also clear that a solution involving light Higgsinos (or a light 
pseudo-scalar) is strongly disfavored by the a s prediction: 

1. The extracted a s line-shape value would decrease (in agreement, however, with low- 
energy extractions). 

2. The predicted a s value would increase due to leading-logarithm [oc — ln(|//|/M^)] and 
possibly non-logarithmic threshold corrections. 
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3. The central value of the fitted m p t ° e [eq. (4b)] would grow to 163 GeV, further increas- 
ing the a s prediction by 0.0003. 

Thus, the Z — >■ bb anomaly, if not resolved, contains strong implications for supersymmetric 
models and could even rule out the simplest and most attractive unification scenarios. 

Lastly, we consider possible high-scale contributions to the correction function. Unlike 
the MSSM case, in which the particles and their mass range are dictated by the model, 
the details of the high-scale corrections are ambiguous. In the minimal SU(5) model [29] 
negative threshold corrections in (7) due to super-heavy color triplet Higgs supermultiplets 
are strongly constrained by the non-observation of proton decay [30], and the GUT-scale 
threshold correction contribution to A„, is typically positive. (This observation, however, 
need not hold in extended models.) Nevertheless, one cannot extract strong constraints on 
the GUT spectrum. Gravitationally induced operators (suppressed by Ma/Mpi anc k ~ 0.001) 
split the Mq gauge couplings (in a correlated manner) and correct the a s prediction in 
proportion to their effective strength, r/, which is a free parameter and can have either sign. 
One has 7 A™° 0.005??. Constraining the NRO corrections to stay perturbative so that 
the calculation is consistent (higher-order terms are negligible in this case) one has \rj\ < 2 
(|r/| ~ 3 is an extreme but still acceptable choice). Thus, NRO's with a non-negligible 
and negative rj could smear light and heavy threshold corrections. Unless rj > and/or 
Msusy *C Mz (which could also imply positive non-logarithmic corrections), no significant 
constraints can be placed on the super-heavy spectrum at present. On the other hand, the 
minimal SU (5) model (where threshold corrections are strongly constrained) would require 
NRO's with j] < if a s (Mz) ^ 0.125. (A similar observation was made recently in Ref. 
[31].) Thus, unification and quantum gravity may be inseparable. 

Regarding the unification scale, corrections that increase the unification scale would 
typically also increase the prediction for a s (Mz) [32], and are thus difficult to construct 
[in particular, for s 2 (Mz) ~ 0.2316]. This is true for contributions to A„, as well as for 
an additional matter family [a s (Mz) — > 0.132] or additional pairs of Higgs doublets [which 
lead to non-perturbative values of a s (Mz)]- This is easily understood if we write a° L as a 
function of the unification scale M and of a(Mz) (see Fig. 4), 



where t = (1/2tt) \n(M/M z ). Naively substituting, e.g., M = M string ~ 5 x 10 17 GeV [33], 
one has a s (Mz) > 0.2. By carefully adjusting operator and super-heavy threshold correction 
contributions to A„ s , one could increase Mq by an order of magnitude while maintaining an 
acceptable prediction for a s [4,8]. However, in general, to rectify the string and unification 
scales (in level-one models) one has to compromise the predictive power of the unification 
scenario [34] so that the correlation between a s and t is modified. 

To conclude, we have shown that typically one expects a large QCD coupling in super- 
symmetric unified models (and even more so when considering a typical MSSM spectrum). 



7 The proportionality factor is calculated here in the SU(5) theory [3], and its normalization is 
different by a factor of four than in [3]. 
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FIG. 4. The Z-pole weak angle and strong coupling are predicted as a function of the unification 
scale M. A given value of s 2 (Mz) corresponds to a fixed choice for M, e.g., s 2 (Mz) = 0.2359 
corresponds to M = 10 16 GeV. MSSM /3-functions are assumed. Two-loop Yukawa corrections are 
taken into account assuming m p t ° le = 160 GeV and tan (3 = 4. (A„, = 0.) s 2 (M z ) = 0.2316±0.0003 
and a s (Mz) = 0.12 ± 0.01 are indicated for comparison. 

This constitutes an interesting signature and has implications for, e.g., Yukawa unification, 
corrections to the unification scale, and the overall sign of the correction function A„ s , and 
is in possible conflict with low-energy data. However, it does not yet allow a significant 
constraint on the super-heavy spectrum because of possible gravitational corrections. We 
also pointed out the interesting role that the Z hadronic width might play in supersymmet- 
ric GUT's, and suggested a simple formula that extends our treatment of mf°' e -dependent 
electroweak corrections to the supersymmetric sector. 
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